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I. INTRODUCTION 



Control system complexity and demands on reliability are 

rising with ever increasing automation. More and more control system 

engineers are looking to solid-state digital techniques for new 

answers to the old problems of reliability, complexity, and 

economy. Digital techniques provide many advantages such as accuracy 
and computational speed as well as flexibility and versatility, 
which analog techniques could not provide. By the use of digital 

techniques in control systems, sensitive, and efficient control 

devices can easily be built. Another significant advantage of 
digital control systems is that the utilizing of digital components 

in a control system allows time sharing of important parts 

of the system which results in economy in the use of equipment. 

The digital elements also require much less power than the analog 
olements. The data signals in digital control systems can also be 
transmitted and received in pulse-code form almost error-free, 

even in transmission through a noisy media. Furthermore, a 

digital computer in control gives another important adventage, it 

makes feasible system compensation by nonlinear programming and 

adaptive or self-optimizing control. 

The digital servo system is becoming widely accepted in 

machine tool control, radar systems, navigational system, and 

many other fields. A digital servomechanism uses numerical 
information to provide a simple but highly accurate method of 

positioning control. Figure 1 shows a block diagram of a typical 
digital servo system. The major components of the basic digital 
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servo system are : 

1. Controlled system 

2. Control elements 

3. Analog-to-Digital converter 
Digital-to-Analog converter 

5. Digital data processing equipment or computer 

The controlled system is normally a linear servomotor but step 
motors have been used for quite some time in both closed loop and 
open loop digital positioning systems. The digital computer is used 

as the central control unit for data processing, system compensa- 
tion, and decision making in some cases. In a simple digital servo 

system, a comparator such as an up-down counter can be used . 

The digital- to-analog and analog-to-digital converters are used to 
interconnect the digital and analog quantities for controlled and 
control elements. One familier type of analog-to-digital converter- 
encoder is the shaft-position encoder disk employing a gray or 

binary code. This digitizer is used to translate the shaft 

position into digital notation by making use of a mechanically 
rotating shaft with coding disk centered on the shaft. Figure 2 

ill ustrates a typical coding disk for a shaft digitizer. The 

brushes contact the disk along the reading line at the point 

shown with a single brush positioned 'in the center 'of each 

concentric band. These concentric rings represent 2^, 2*, 2^,...,2 n , 

and each ring . ' is ' constructed of several segments made of 

either conducting material (the darkened areas) or insulating 

matorial ( the unshaded areas). A signal is connected to the 
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conductor, and if a given brush makes contact with a segment of 

conducting material, a " 1 » signal will result, but if the brush 
is over the insulating area, the output from the brush will be 

a " 0 The resolution of this encoder is determined by the 

number of channels. For a 4-channel coding disk, a resolution of 
one part in sixteen is achieved. Two types of encoder disk, the 

gray and binary code disks are shown in Figure 3. The coded disk 
my be made of transparent and opaque materials arranged in coded 
patterns. The photocells and sources are used to read the angular 

position in a digital code. 

In operation, the analog-to-digital converter converts the 

analog output of the controlled system into a numerical quantity, 
which is compared with a numerical order, and the difference is 
converted by the digital-to-analog converter into an analog signal 

capable of driving a servomotor until the difference is eliminated. 

A. OBJECTIVES 

The methods for analysis and synthesis of the analog system 

such as Bode plots. Root Locus, Parameter plane plots, etc. have 
been established for a long time. These methods work very well, 
however, for the discrete or digital system, there is no such 
method that is as good as those for the analog system. Even 

the well known z-transform method can become very complicated if 

the order of the systems analyzed are higher than third or fourth 
order. 

The objective of this paper • was to investigate a typical 
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digital positioning servo system in order to obtain some good 

methods for analysis and synthesis. In the proposed system, light 

travels from a light source through equally spaced openings in the 

code disk to activate the sensor. The sensor, in turn, generates 
feedback pulses at equilibrium angular position. The sampling 

period of this system is not constant, depending on the velocity 
of the motor, the z-transform method of analysis and design cannot 
be used here. Another objective was to find a suitable mathematic 
representation of such a system. 

All of the studies had been done by simulation using Digital 
Simulation language DSL/360 on the IBM 3&0 computer system, this 

language made the simulation quite easy, however, care must be 

taken in order to obtain accurate integration in the system with 

high velocity. 

In Chapter II, the system identification and method of 

simulation are given. Chapter III presents the discussion of the 
results obtained from the simulation of second and third order 

servo systems. This includes the study of the stability limit, 
percent overshoot, and the effects' of variation oof the 

servoamplifier gain. This chapter also shows the results obtained 
from investigating the behavior of the system with ramp input. 
The Bode plot and parameter plane plot were used to compensate 

the analog system in Chapter IV. The results were then used in 
simulating the digital system, since there was no direct compensation 
method. 

The results were obtained both in graphs and printed output, 
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however, only graphs are included. Listing of all computer programs 
are presented after the conclusion and recommendation in Chapter V. 

In each study, the digital system with several numbers of 
opening slots and the continuous counterpart were simulated. The 

studies of the digital systems were based on the results of studies 
of the continuous system. 
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II. THE DIGITAL POSITIONING SERVO 



A. PHYSICAL PROPERTIES 

The digital positioning servo system under consideration consists 
of a digital comparator, a digital -to-analog converter, a controlled 
system, and an analog -to-digital converter. The analog -to-digital 
converter consists of a low inertia encoder disk mounted tightly 
on the servomotor shaft, two photoelectric sensor units, and a 
direction-sensing logic circuit. The encoder disk has uniformly 
spaced slots around its circumference. Each photoelectric sensor 
unit consists of a light source and a photodiode. Figure 4 

illustrates the encoder disk and the sensor units and Figure 5 
illustrates the direction-sensing circuit. 

Each time an opening passes through the beam from the light 
source, a pulse is generated by the sensing unit, and it will 
eventually be fed to the up-down counter which acts as the 

digital comparator. The second sensor unit is for the purpose of 
direction sensing. When the servomotor turns in forward direction, 
the feedback pulse becomes a down-enable signal, when the servomotor 
turns backward, the feedback pulse becomes an up-enable signal. 

The two sensors are adjusted so that their respective outputs, X 
and Y, overlap as shown in Figure 6. Figure 6-A and 6-B show 
the inputs and outputs of the direction-sensing circuit when the 
servomotor turns forward and backward, respectively. Figure 7 shows 
a typical digital-to-analog converter. 

For an input as a step function, an appropriate binary number 
representing the size of the step function is fed into the up-down 
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counter. The output of the up-down counter is, in turn, converted 
into an analog signal and is sent to the servoamplifier to 

drive the servomotor. The sensing device will measure the shaft 
position, encode, and send this information to the up-down counter. 
The positioning process will be finished when the binary number 
in the up-down counter becomes zero. 

B. COMPUTER SIMULATION 

The DSL/360 language was specially designed for use in 

simulation. D5L/360 provides . many advantages over other languages 
and subroutines in existence, especially in simplicity and time 
saving. DSL/360 not only provides these advantages in simulating 

the linear physical system, but in simulation of the non-linear 
system also, DSL/360 provides a basic set of function blocks 
from which a physical system may be modeled such as integrators, 
quantizer, pulse generators, function generators, and so on. In 
addition, FORTRAN Library functions from the Scientific Subroutine 
Package may be utilized. The user may also provide his own 
subroutines if necessary. 

In programming, first, each part of the digital positioning 
servo system was separated into single function blocks as shown 



in Figure 8. 


Next, 


each 


single function block 


was 


replaced by a 


related DSL 


f vine t ion 


or 


standard FORTRAN Library 


function as 


appropriate. 


Finally, 


all 


the DSL and FORTRAN 


Library function were 



put together- in the DERIVATIVE REGION of the program. One advantage 
obtained in using the DERIVATIVE REGION is that the functions are 
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not required to be in any specific order. 

In the early parts of simulation. The servomotor was assumed 
to be a second order motor with the Transfer function 



K v 

G (s) = 

m s(s + 1) 



Later, the third order servomotor with the transfer function 



G (s) 

s (s + 1) (0. 2s + 1) 

was used. The FORTRAN Library function 

Y = AINT(X) 

was used to modelled the quantizer. This function was also used 

to quantize the ramp input function. The gain of the digital- 

to-analog converter was designed to be equal to l/K to make the 

c 

system a unity feedback system. 

The fourth order Runge-Kutta with fixed integration interval 
(INTEG RK3FX) Integration routine was used in all simulations to 

assure the accuracy of the results. It is very important to 

choose a correct integration interval (BELT), since too large DELT 
will cause inaccurate integration and too small DELT will result 
in excessive computing time. The DELT used in each simulation 
is approximately five to ten percent of the time interval between 

two consecutive pulses generated by the digitizer at its maximum 
angular velocity in that specific simulation. 

In multiple runs, the TERMINAL region was used to terminate 
the previous run, to provide the new sots of parameter values. 
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and to terminate the program after the output data of the last 
run had been plotted. The FORTRAN subroutine DRV/G was used in 

the SAMPLE region to collect points for each sample interval 

(DELS) to be plotted by the CALCOKP XY Plotter. 

In most cases. The digital system with 10, 45, 90, and 180 

openings and the continuous system were simulated and the results 
were compared. A system with one opening was used in stability 
limit studies since it was the worst case. To compare the outputs 

of the digital system with output of the continuous system, the 

inputs of all systems must be equal. If a step input to the 

continuous system is 1.25664, inputs to the digital system with 
10, 45, 90, and 180 openings must be 2, 9.-18, and 36 respective- 

ly. 

These input values can be readily calculated by the method 

shown in example below. 

EXAMPLE 

To calculate comparable step inputs to the continuous and 
digital system with 10 openings. 

From Figure 8, 

^in( continuous) ” ^ ^ 

R - IGBIT 

V in(digital) k 

c 

2 (R n - IGBIT) 

OP 



for zero initial conditions, 

^in ( continuous ) 
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2 R n 

V in(digital) = 

to have equal inputs to servomotors, requires 

2 R 

R = £ 

OP 

by arbitarily choosing R^ = 2 

R = 1.25664 

The values of inputs to digital system with other numbers of 
opening were also based on this selection. 
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III. DISCUSSION 



The first part of this chapter contains a discussion of the 

response of the second order servo system to a step function 

input. The third order system with both step and ramp function 

input is discussed in the later part of the chapter. 

A. STEP RESPONSE OF THE SECOND ORDER SYSTEM 

The main concern in the response of the second order servo 

system was its stability limit. The forward gain of the system 

had been increased until the stability limit was obtained. 

The model of the second order continuous system used in this 

simulation is always stable regardless of the size of the gain. 

The digital system was also stable at all value of gain regardless 

of the number of openings in the position encoder. The gain 

7 

of the system was increased to 10 to make sure that the system 

would always be stable. Figure 9» 10, 11, 12, 13, and 14 show the 
step response and plots of the binary numbers in the up-down 
counter versus time for the digital system with number of openings 

equal to 1, 10, 45, 90, and 180, and for the continuous system, res- 
pectively. Figure 15, 16, 17, 18, 19, and 20 show their respective 
phase plane plots. The plots of the output of the up-down counter 
were not scaled correctly, however, they could be read oasily. 
Each quantizied step of the N^ plot represents a one. 

The difference that can readily be seen in the responses of 
the digital and continuous systems is the magnitude of both the 
angular displacement and velocity. The difference is very pronounced 
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when the number of openings is small. The digital system with 
fewer openings has longer time delay between two consecutive 

feedback pulses. This time delay has allowed the motor to move 
farther with higher velocity, similar to the response of the 
system without a feedback loop, before the excitation signal is 
reduced by the next feedback pulse. The system with a large 
number of openings has a very high sampling rate, therefore this 

system behavior is much closer to that of the continuous, system, 

a code disk with 20000 openings has been developed. This disk 

gives the system almost the same response as that of a continuous 
system. 

The phase plane plots of the digital system are not as 
smooth as those of the continuous system. They actually consist 
of many short sections. Each portion represents the delay 

interval between two consecutive pulses. This can easily be observed 
in the phase plane plot of the system with lower number of 
openings such as shown in Figure 16. 

Figure 21 contains the comparison of the angular displacement 

of the analog and digital system with 10 and 45 openings. This 
plot clearly shows that the natural frequency of this second 

order digital system is smaller than that of the analog system. 

It indicates that the root locus of the digital system would 

cross the CO -axis closer to the origin than the root locus of 
the continuous system. This distance becomes smaller if the number 
of openings is decreased. 
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The plots of the step response of the digital system also 
consist of many short sections. It is very easy to observe this 
phenomenon in the response of the system with very low gain 

such as those shown in Figure 22 and 23. Figure 24 demonstrates 
how this phenomenon occurs. At the beginning of the simulation, 
the system behaves like a system without a feedback loop. The 

rate of increase of velocity is very high at first, however, it 
becomes lower with time. The displacement plot also becomes more 
linear with time. When the displacement' is large enough, a feedback 
pulse is generated, the output of the up-down counter is changed 
by one, the velocity and the rate of change of the displacement 
are changed to follow the nGW input voltage to the servoamplifier. 
This process occurs repeatly until the output of the up-down 
counter is eliminated, and the system eventually comes to rest. 

Figure 22 and 23 also illustrate that the digital system behaves 
like the continuous system in another aspect. When the gain is 
reduced, the peak overshoot decreases. This peak overshoot can be 
eliminated by making the gain small enough. It is also observed 
that, for all cases studied, the steady state displacement of the 

digital system is always larger than that of the continuous 

system with comparable size of step input. The difference, however, 
has never exceeded the angular displacement between two adjacent 
openings. The system with more openings in the encoder disk gives 

smaller error than the system with fewer openings. The maximum 

error for the system with 10, 45, 90, and 180 openings are 

less than 0.6283, 0 . 1396 , 0.0698, and 0,0349 radians respectively. 
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B. THIRD ORDER SYSTEM 

A third order system was studied in greater detail than 
the second order system. In addition to what had been studied 
from the second order system characteristics, the variation of 
the percent overshoot was investigated in this part. 

1. Stability limit 

A third order continuous system with the transfer function 
shown in Chapter II reaches the stability limit when the 
servoamplifier gain, Ky is equal to 6. The digital systems, 
however, require larger gains to reach their stability limits. 

Figure 25, 26, 27, 28, and 29 illustrate the phase plane plots of the 
third order digital system with the gain equal to 6 and the 

number of openings equal to 1, 10, 45, 90, and 180 respectively. 

Figure 30 shows the phase plane plot of the continuous system 

with the same gain. The phase plane plots of the third order 
digital system converge into a limit cycle and can clearly be 

observed in the case the system has 1 or 10 openings. By 
inspecting the characteristic of the phase plane plots, especially 
in the case where the number of openings is small, it is 

obvious that the third order digital system requires higher gain 

to reach the stability limit than the continuous system. The 
time delay also causes tho displacement, velocity, and limit 
cycle size of the digital system to be larger than those of tho 
continuous system. 

Since there had not been any known method that can bo 

used to find the gain at stability limit of this digital servo 
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system, a trill and error method was used to find this gain, the 
gain of the third order digital system with 10 openings at its 
stability limit was found to be approximately 6.25. This gain 
would be reduced toward 6 if the number of openings is increased. 
The size of the limit cycle is also decreased with the increase 
of the number of openings. 

The results from simulating the third order digital servo 

system with number of opening equal to 10 as illustrated in 

Figure 31 * 26, 32, and 33 also show that the size of the limit cycle 
increases with the gain. The increase of the gain also causes 
the system to reach this limit cycle faster. Upon making the 

gain greater than the stability limit gain, The phase plane plot 
diverges for a period of time before reaching another limit cycle 
as shown in Figure 34-. For higher gain, The phase plane plot 
diverges into a limit cycle and stays in the limit cycle for a 
period of time then diverges to infinity as shown in Figure 35. 
With the gain a little higher as shown in Figure 36, the phase 

plane plot diverges to infinity but has a tendency to go into 

a limit cycle at first. When the gain is much higher, the 

system just becomes unstable as shown in Figure 37. 

2. Percent overshoot 

For a continuous and stable servo system, the overshoots 
of the step response always decrease with time, and the plot of 
the log of the percent overshoot versus time is a straight 

line. The step response and the plot of log of percent overshoot 
of the continuous system are shown in Figure 33 and 43 respectively. 
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To investigate the behavior of the overshoot of the digital 
positoning servo system, a gain of 2.5 was used. With this gain, 

£ and oj n of the continuous system are equal to 0.166 and 1.54 

respectively. Figure 39, 40, 41, and 42 illustrate the step response 
and output of the up-down counter of the third order digital servo 

system with 180, 90, 45, and 10 openings respectively. The step response 
of the continuous system shows 6 prominent overshoots, the responses 

of the digital system with 180, 90, 45, and 10 openings show 5, ,4, 
3, and 2 overshoots respectively. The time at which each of the 

corresponding peaks of all systems occur is the same. Figure 43 

illustrates the plots of the percent overshoot of the responses 
of the continuous and digital system with 180, 90, 45, and 10 

openings versus time. The numerical values of percent overshoots 

are included in table I, The inputs to all systems were adjusted 
so that they would be equal and the percent overshoot could be 
calculated from the same basis, the steady state output of the 
continuous system. 

The overshoots of all digital systems are higher than the 
corresponding peaks of the continuous system. The difference, 

however, does not exceed the angle of the two adjacent openings. 
The digital system with smaller number of openings has larger 

percent overshoot than tho system with larger number of openings. 

The steady state response of the digital system is also larger 

than that of the continuous system. The difference is also less 

than the angular displacement between two adjacent openings. 

At the beginning, the plot of tho percent overshoot of the 
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digital system tends to be parallel to that of the continuous 
system. The plot of the percent overshoot then diverges to the 
smallest peak which is slightly larger than the angular displacement 
between two adjacent openings. Since the system with fewer 
openings has larger displacement between the openings, its response 
requires less time to reach the last overshoot. The number of its 
overshoots is also less than that of the system with more openings. 

The responses of all the digital systems that have less 
than 10 openings were also studied. All of the responses have 
two overshoots and the characteristics are similar to those discussed 
above. However, only the response of the system with one opening, 
which is the extreme case, is illustrated in Figure 44 . 

3. Response to ramp inout 

The ramp input function, Bt, was also applied to the third 
order digital positioning servo system in order to investigate its 
stability. The servoamplifior gain used in this study was 2.5 
and the slope of the ramp was varied from 0.5 to .5. . Practically, 
a method of applying a ramp input to the digital servo system 
is to apply a series of impulses into the up-down counter. The 
slope of the ramp can bo controlled by varying the period of the 
impulse. A digitized ramp function was used ' in this simulation 
instead of the impulses. The study of the response of the digital 
servo system with ramp input was done by simulating the system 
with 10 openings only. It had been shown from the previous studios 
that the digital system with more openings behaves closer to the 
continuous system than the systom with fewer openings. 
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Figure 45 and 46 illustrate the response of the digital system 
with 10 openings and the slope of the ramp input equal to 0.5, 

1.0, 1*5, 2.0, 2.5, 3.0, 4.0, and 5.0. When the slope is small, 

for example B = 0.5, the rate of the input pulse applied to the 

up-down counter is slow and causes the velocity of the shaft and 

the feedback sampling rate to be too slow and nonuniform. The 
nonuniform and slow rate of output of the up-down counter, in 
turn, cause the shaft to turn abruptly and the system becomes 
unstable. The longer time delay is also an important cause for 

the instability of the system. The plot of velocity and output 
of the up-down counter for the digital system with the slope of 

the ramp input to the system equal to 0.5, 1.0, 1.5, 2.0, 2.5, 

3.0, 4.0, and 5.0 are shown in Figure 4 7„ 48, 49, 5Q, 51, 5?, 53, 

and 54 respectively. When the slope of the ramp input is increased, 

the input pulses are fed to the up-down counter at higher rate, 

the servo shaft turns faster, the feedback sampling rate becomes 
higher, the output of the up-down counter becomes more linear and 
the velocity of the shaft becomes more stable. Unlike the continuous 
system, the velocity of the digital system will not be constant 
no matter how large a slope is used for the ramp input. It rather 

varies within a small limit. This limit cycle is caused by the 

time delay between pulses. The limit cycle will be smaller if 

the time delay is shorter. 

Another interesting result observed ' is that the rate of 

decrease of the magnitude of the limit cycle with the increase 
of the slope of the ramp is nonuniform. At first the magnitude 
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of the limit cycle becomes very small when the slope is increased 
to 2.5. However, when the slope is increased to 3 and 4, the 

magnitudes of the limit cycles are some what larger. The limit 
cycle is reduced to a smaller size again when the slope is 

increased to 5. 
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IV. SYSTEM COMPENSATION 



Since there had not been any known method that can be used 

to compensate this digital positioning servo system directly, the 
author chose to compensate the comparable continuous system with 
known methods and simulated the compensated . system as both a 

continuous and a digital system. The results from the simulation 

were then compared and studied. The continuous and the digital 

system with 10, 45, and 180 openings were simulated and their 

step responses were plotted on the same graph to make the 
comparison more obvious. The simulation of the digital system with 

90 openings was omitted to prevent the graph from being overcrowded. 
The system compensated was the third order servo with servoamplifier 

gain equal to 6. 

A. LAG COMPENSATION 

Since the programs and computer facility were available, the 

parameter plane plot method was used in lag-compensating of the 
system. This method is very efficient and flexible and the computer 

make it very easy to obtain the plot. The author selected the 

PARAK K program to use in this compensation. The complete details 
of this program are shown in the appendix. 

In order to obtain the parameter plane curve, the characteristic 
equation of the system and other requirements such as the value 
of £ , ton , Co , or d needed to be plotted and the scales of the 
plot must be provided as the input data to the PARAK K program. 

Figure 55 illustrates the block diagram of the third order 
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continuous system with the lag compensation section. The lag 
compensation block has the transfer function 



G 




* D 

(T p s + i) 



The transfer function of the third order system with 
compensation becomes 



G (s) 

m 



(T z s +D 6 

■ 0 ■ "" — ■ ■ ' ■ u .1 

(TpS + 1) s(s + 1) (0.2s + 1) 



the lag 



or 



G(s ) --- £ < 1 1 z) 30 

z (s + p) s( s + 1) (s + 5) 

The characteristic equation of this system is 

s 4 + (6 + p)s 3 + (5 + 6p) s 2 + (5p + 30| )s + 30 P = o 



let oc = p 




the characteristic equation becomes 



s 4 + (6 +a)s 3 + (5 + 6a)s 2 + (5 •f-30f3)s+30a:=0 

The parameter plane plot of this characteristic equation is shown 
in Figure 56 . oc and 0 are the unknown pair of variables which 
their values can be directly selected from the suitable value of 
5 and GJn. The values of pole and zero can also be readily 



obtained from the values 


of oc and /3 


. Three pairs of 


poles 


and 


zeroes were chosen and 


simulated. 


The 


step responses 


of 


the 


compensated continuous 


and digital 


system 


with 10 , 45 , 


and 


180 
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openings are illustrated in Figure 57, 58, and 59 . The values of 
pole, zero, C, GJn , and predicted and resulted values of M t , 

pt p 

and t of the compensated system are included in table II. 

•The simulated results show that the characteristics of the 
digital system are similar to those of the • continuous system. 

The results also agree with those studied in the previous chapter. 
The magnitude of the peak overshoot and steady state response of 
the digital system are always larger than those of continuous 

system. All the differences are also within the limit of the 
angular displacement between two adjacent openings. 

Care must be taken in selecting the value of the zero for 
the lag compensation section. If this zero was too close to the 
origin, the real root would dominate the system instead of the 
selected pair of imaginary roots. If the zero was too far 

from the origin, the value of 5 would be too small and cause 

the system to be very slow. 

B. LEAD COKPENSATION 

The same characteristic equation was used to obtain the parameter 
plane plot for lead compensation. However, the scales, the starting value 
of Cdfi , the required Values of ? and U>n curves, and etc. must 

be reselected in order to move the curves into required region 
for the lead compensation. The block diagram of the system with 
lead compensation section is the same as that of the lag 
compensated system. Its parameter plane plot is shown in Figure 60. 

The pole and zero of the lead compensation section were kept 
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within one decade apart to prevent the severe reduction of the 
forward transfer gain of the section. This imposes another limitation 
into the selection of the value of (3. The value of £ should 
not be larger than 10. Three pairs of poles and zeroes were 

selected and simulated. Figure 6l, 62, and 63 illustrate the step 
response of the compensated system in both continuous and digital 
versions. Table III consists of the values of pole, zero, 5 , cj n . 



predicted 


and 


resulted values of t , 

P 


V 


and t of the 
s 


compensated 


system. 


Since 


the natural frequency 


of 


the lead 


compensated 



system is larger than that of the lag compensated system it 

results in smaller t and t . The other results of this study 

p s 

are similar to those described in the lag compensation study above. 

C. LAG -LEAD COMPENSATION 

The block diagram of the third order continuous system with 

lag and lead compensation sections is shown in Figure 6k. There 
are four variables in this system which make it more difficult 
to use the parameter plane plot method in compensation. The Bode 
plot method was used in this compensation. This method is rather 

simple but requires some degree of neatness . Figure 65 and 66 

illustrate the Bode plots of two different lag-lead compensations-. 

The limitation in selecting the values of the zero of the lag 
section and the spread between the pole and the zero of the lead 

section are also imposed here. Both of the compensated systems 

were also simulated and studied. The results of the simulation 
of the lag-lead compensated systems are shown in Figure 67 and 68„ 
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Table IV consists of the values of pole, zero, f , cj n , predicted 
and resulted value of t , M ., and t . The results of the lag- 

P pT/ S 

lead compensation study are also similar to those of lag and lead 
compensation studies, 

D. COMPENSATION DISCUSSION 

The results from all compensation studies as well as those from 
the previous chapter show that the steady state response of the 

digital system is always within the angular displacement between 
two adjacent openings larger than that of the comparable analog 

system. However, this difference is not constant, it changes randomly 
within the limit. The exact position of the steady state response 

of the system can not be determined. It seems that the best method 
to compensate the digital system is to compensate its analog 

counterpart. 

The maximum error of the system under investication is found 
to be 

1.296X10 seconds 
OP 



In comparison, the analog system has the accuracy of approximately 
2 minutes of arc and the digital system with gray or binary 

shaft-position encoder has the accuracy of 



1.296X10 



6 



2 



n 



seconds 



where n is the number of the conducting channel (see ref. 5). In 
order to have the same accuracy as of the analog system, this 
digital system must have approximately 10000 openings in the encoder 
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disk. The digital system with gray or binary code disk requires 
13 channels in the code disk to acquire comparable accuracy. It 
is trivial that this digital system has an economic advantage over 
the gray or binary code disk system. This system requires only 

two sensing units for any number of openings. The gray or binary 

shaft-position encoder unit requires one sensing device for each 

channel in the code disk and the sensing circuit has to be more 
complicated. 

Another adventage is that all the studies show that the 

response of this digital system is always larger than the response 
of the continuous system. This observation could be used to reduced 
the steady state error by making the size of the step input to 

be less than the required input to the system. The difference 

should be as close to a half of the maximum error as possible. 
The difference of exactly a half of the angular displacement 

between two adjacent openings would reduced the steady state error 
to within one half of the maximum error. 
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CONCUSSION 



The characteristics of the digital positioning system under 

investigation can be very similar to the analog system if the 

feedback sampling rate is high enough. This system has some 
characteristics which are superior to the conventional' analog and 
sampled-data system, especially in accuracy and economy. 

Compensation of this system can be accomplished by using the 
results obtained from compensating its analog counterpart by 

conventional methods. The accuracy can be readily improved by 
modifying the encoder disk. 

From the discussion above, it can be seen that there are 

many areas left to be studied. The author believes that 
further study of the mathematical analysis and of the direct method 
of compensation are necessary. 
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Figure 1. Block diagram for a typical digital control system. 




Figure 2. A typical shaft position-to-digital converter 
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binary code gray code 

Figure 3. Two typical code disks. 




Figure 4. Code disk and photoelectric sensor units. 







Figure 5. Direction-sensing circuit. 
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Figure 6. Input and output signals of the direction-sensing circuit 
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Figure 7. Digital-to-analog converter 




a) continuous system 




b) digital system 



Figure 8. Positioning servo system in 
single function block diagram 
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Figure 9. Second order digital servo, 1^=10 ,,R n =l, OP=l 
displacement and output of U/D counter vs time 
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Figure 10. Second order digital servo, K v =10 , R n = 2 * • 0P=i0 

displacement and output of U/D counter vs time 
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displacement and output of U/D counter vs time 
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displacement and output of U/D counter vs time 
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displacement and output of U/D counter vs time 
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displacement vs time 



Velocity, radians/second 
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Figure 15. Second order digital servo, K v =10^, R n =l, 0P=1 

phase plane plot 
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Velocity, radians/second 




phase plane plot 
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Figure 17. Second order digital servo, K =10^, R =9, 0P=^5 

I V n 

phase plane plot 
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Velocity, radians/second 




Figure 18. Second order digital servo, K V =10^, R n =18, 0P=90 

phase plane plot 
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Figure 19. Second order digital servo, K v =10 , 

phase plane plot 



R = 36 , CP=180 
n 
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Velocity, radians/second 




n 

Figure 20. Second order analog servo, K v =10', R=l. 25664 

phaso plane plot 
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Displacement, radians 




-^Figure 22. Step response of second order digital servo 

K si .5, l.o, 0.5, 0.36, 0 . 25 , 0.1, and 0 . 05 , 

0P=10, and R n =2 
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Displacement, radians 




•Figure 23. Step response of second order digital servo 
K' v =1.5» 1.0, 0.5, 0.36, 0.25, 0.1, and 0.05, 
0P=45, and R n =2 
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Velocity (2), radians/s ocond 




Time, seconds 



Figure 24. 



Second order digital servo, 1(^=0. 05, OP=10, and R^=2 

Up-down counter output, velocity, and displacement vs time 
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Figure 25. Third order digital servo, K v =6, R^=l, 0P=1 

phase plane plot 
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Velocity, radians/second 




Figure 26. Third order digital servo, K =6, R =2, 0P=10 

phase plane plot 
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phase plane plot 
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Figure 28. Third order digital servo, K y =6, R^=18, OP=90 

phase plane plot 
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Figure 29. 
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Displacement, * t Kio 1 ~ x radians 

Third order digital servo. K =6. R = 36 . 

v n y 9 

phase plane plot 
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Velocity, radians /second 




Figure 30. Third order analog servo, K^=6, R=l. 25664- 

phase plane plot 
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^.oc ij 2o.:»: 

Displacement, :xio ? radians 

Figure 31. Third order digital servo, K v =5.7, R n =2, OP=10 

phase plane plot 
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Figure 32. Third order digital servo, K v =6.1, ^ n =2, 0P=10 

phase plane plot 
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Figure 33. 



Third 



order digital servo, K v =6.25, R n =2, 



phase plane plot 



0P=10 
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Figure 3^. Third order digital servo, K v =6.5> R n =2, OP=lO 

phase plane plot 
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Figure 35 . 



Third order digital servo, K v =7.1, R n =2, OP=lO 
phase plane plot 
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Displacement, radians 

Figure 36 . Third order digital servo, K v =?.2, R n =2, OP=10 

phaso plane plot 



7 ? 



Velocity, radians/second 




Displacement, radians 



1 

.V. 



Figure 37. Third order digital sorvo, K v =7.5. R n =2, 0P=10 

phase plane plot 
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Displacement, radians 




Time, seconds 



Figure 38. Third order analog servo, K^=2.5, R=l. 25664 

displacement versus time 
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Displacement (1), radians 
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Figure 39. Third order digital servo, K =2.5, R =36, 0P=180 

v n 

displacement and up-dovm counter output vs time 
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Up-down counter output (2) 



Displacement (1), radians 
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Figure 40. Third order digital servo, K v =2.5, R.=18, OP=90 
displacement and up-down counter output vs time 
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Up-down counter output (2) 



Displacement (1), radians 




Time, seconds 



Figuro 41, Third order digital servo, K v =2.5, R^=9. 0P=45 
displacement and up-down counter output vs time 
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Up-down counter output (2) 



Displacement (1), radians 




Time, seconds 



Figure 42. Third order digital servo, K v =2.5, R n =2, CP=10 
displacement and up-down counter output vs time 
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Up-down counter output (2) 



Percent overshoot 




Figure 43 . Third order servo system, K =2.5 
OP=10, 45, 90, 180, and analog 
percent overshoot vorsus time 
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Displacement, radians 




Time, seconds 



Figure W*. Third order digital servo, K v =2.5» F =2, 0P=1 

displacement versus time 
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Displacement, radians 




Time, seconds 

Figure 45. Third order digital servo, displacement vs time 
ramp input, B=0.5, 1.0, 1.5, and 2.0, OP=lO, K y =2.5 
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Displacement, radians 




Figure 46. Third order digital servo, displacement vs time 
ramp input, 3=2.5, 3.0, 4.0, and 5.0,,0P=10, K v =2.5 
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47. Third order digital servo, ramp input, B=0.5, K y =2.5, OP-10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians/second 




Figure 48. Third order digital servo, ramp input, B=1.0, K^=2.5, 0P=:10 
velocity „ and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians/second 







Figure ^9. Third order digital servo, ramp input, B=1.5» K v =2.5. OP=10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians/second 




Figure 50. Third order digital servo, ramp input, B=2.0, K v =2.5, 0P=10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians /second 




Figure 51. Third order digital servo, ramp input, . B=2, 5, K v =2,5, 0P=10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians /second 




Figure 52. Third order digital servo, ramp input, B=3.0, K v =2.5, 0P=10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians/second 




Time, radians 

Figure 53. Third order digitalservo, ramp input, B=4.0, K v =-2.5» OP=10 
velocity and up-down counter output versus time 
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Up-down counter output (2) 



Velocity (1), radians/second 




Figure 54. Third order digital servo, ramp input, B=5.0, K y =2 . 5 , OP=10 
velocity and up-down counter output versus time 
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with a lag compensation section 



C * \ o 




Figure 5^. Parameter plane plot for lag compensation 
of the third .order positioning servo 
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Figure 5?. Lag compensated third order positioning servo 
P=0.013, z*=0. 18, 0P=10, 45, 180, and analog 
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Displacement, radians 
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Figure 58 , Lag compensated third order positioning servo 
p =0.017, z=0.19, 0P=10, 45, 180, and analog 
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Figure 59. Lag compensated third order positioning servo 
p^O.Ol, z=0.15, 0P=10, ^5. 130, and analog 
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Figure 6 q. Parameter plane plot for lead compensation 
of the third order positioning servo 



101 



i>isplacement, radians 






8 

ki 
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Figure 6l. Lead compensated third order positioning servo 
P=14.5, z~1.85i 0P=10, 45, 180, and analog 
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Displacement, radians 




Time, seconds 



Figure 62. Lead compensated third order positioning servo 
P=15.0, z=l,50, 0P=10, 45, 180, and analog 
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Displacement, radians 




Time, seconds 



Figuro 63 . Lead compensated third order positioning servo 
p=l9.0, z=1.90, OF=lO, 45, 180, and analog 
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Figure 64. 31ock diagram of a third order servo 
with a lag-lead compensation section 



qp ‘up t?q 
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Figure 65. Lag-lead compensated third order servi 

Plag^* 25, Z lag =0,5, p lead =20, Z lead =: 
Bode • plot 
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Figure 66. Lag-lead compensated third order servo 

Plag^* 2, Z lag =0# ' 5, p lead =12,5, Z lead“ 1,25 
Bodo plot 



Displacement, radians 



r f 
rj 




~l ' I 1 1 1 

■•..OG 1.03 12. 03 :<v03 . 2. 

Time, seconds 



Figure 67 . Lag-lead compensated third order servo 

p, =0.25. z. =0.5, p, .=20 , z. ,=2 
r lag lag 1 lead ’ load 

OP = 10, 45, 180, ,and continuous 



108 



Displacement, radians 




Time, seconds 



Figure 68. 



Lag-lead compensated third order servo 

Plag"" 0,2, Z lag =0,5, p load =12 ‘ 5, Z load =1 * 25 
OP = 10, 45, 180, and continuous 
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APPENDIX 



This appendix presents four examples of the programs used to 
simulate the system under investigation and a program used to plot 
the parameter plane curves. 
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TITLE SECOND ORDER ANALOG POSITIONING SERVO 

* STEP RESPONSE AND PHASE PUNS PLOT 

* 



* 

I NT EG RKSFX 
INTGER NT LOT, NC UR 

CONST NPL0T=2 , NCUR=2,KV=1,0E?,R=1. 25664 
INITIAL 
DERIVATIVE 
VT-R-C 
CDD=VT*KC 

CD0T=REALFL(0.0,1.0,CDD) 

C=INTGRL ( 0 . 0 , C DOT ) 

CONTRL FINTIM=0. 01 , DELT=0 . 00000 1 , DELS=0. 000002 

PRINT 0.0001 ,VT,CDOT,C 

SAMPLE 

CALL DRWG(1,1TIME,C) 

call drv;g(2,i,c,cdot) 

TERMINAL 

CALL ENDRW ( NPLOT) 

END 

STOP 



HI 



TITLE SECOND ORDER DIGITAL POSITIONING SERVO 

* STEP RESPONSE AND OUTPUT OF THE UP-DOWN COUNTER 

* 

* A IS A FACTOR USED TO ADJUST THE SIZE OF THE CURVE 

* 

* MULTIPLE RUNS 

* 

I NT EG RK3FX 
INTGER NFLOT , NCUR 

CONST NPLOTsl , NCUR=5 , OP=l , RN=1 , KV=1 . OE? , A=20 . 0 
INITIAL 

KC=OP/6. 2831 853 
DERIVATIVE 

BIT=C*KC 

IGBIT=AINT(3IT) 

NT=RN-IGBIT 

VT=NT/KC 

CDD=VT*KV 

CDOT=RSALPL( 0 . 0 , 1 . 0 , CDD) 

C=INTGRL ( 0 . 0 , C DOT ) 

NO=NT*A 

CONTRL FINTIM=0 . 01 , DKLT=0 . 000001 , DELS=0. 00002 

PRINT 0.0001 ,VT,CDOT,C, NT 

SAMPLE 

CALL DRWG ( NPLOT , 1 , TIME , NO) 

CALL DRWG ( NPLOT , 2 , TIMS , C ) 

TERMINAL 

IF(NCUR.SQ.l) GO TO 6 
NPL0T=NPLCT+1 
NCUR=NCUR-1 
IF(NCUR.E3.4) GO TO 1 
IF(NCUR.22.3) GO TO 2 
IF(lvCUR.E3.2) go to 3 
IF(NCUR. EQ. 1) GO TO 4 

1 0P=10.0 
RN=2.0 

A=1.5 

GO TO 5 

2 0P= 45.0 
RN=9.0 
A=0.3 
GO TO 5 

3 OP=90.O 
RN=18.0 
A=0. 15 
GO TO 5 

4 0P=180.0 
RN- 36 .O 
A=0.075 

5 CALL RERUN 
RETURN 
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6 CALL ENDRVf(NPLOT) 



END 

STOP 



. 1.13 



TITLE THIRD ORDER DIGITAL SERVO 

* VELOCITY AND OUTPUT OF UP-DOWN COUNTER 

* VARIATION OF SLOPE OF RAMP FUNCTION 

* 

* MULTIPLE RUNS 

* 

I NT SO SK3FX 
INTGER NPLOT, NCUR 

CONST NTLOT=l , NCUR=4 , KV=2 . 5 , OP=lO . 0 , 3=5 . 0 
INITIAL 

kc=op/6 . 2831853 

DERIVATIVE 

RA!-j=3*RAi-P (0.0) 

RN=AI’JT(RAM) 

BIT=C*KC 

IGBIT=AINT(BIT) 

NT=RN-IGBIT 

VT=NT/KC 

CD=VT*KV 

CDD=REALPL( 0 . 0 , 0 . 2 , CD) 

C DOT=REALPL ( 0 . 0,1.0, CDD) 

C=INTGRL ( 0 . 0 , C DOT) 

NO=NT 

CONTRL FINTIN=30 . 0 , DELT=0 . 004, DELS=0 . 03 

PRINT 0. 2, VT,CDOT,C, NT 

SAMPLE 

CALL DRUG ( NPLOT , 1 , TIKE, NO) 

CALL DRUG ( NPLOT , 2, TIPS, C DOT) 

TERMINAL 

IF(NCUR.EQ.l) GO TO 6 
NCUR=NCUR-1 
NPL0T=HPL0T+1 
IF(NCUR.EQ.3) GO TO 1 
IF(NCUR. S3. 2) GO TO 2 
IF(NCUR.EQ.l) GO TO 3 

1 B=4, 0 
GO TO 5 

2 B=3.0 
GO TO 5 

3 B=2 . 5 

5 CALL RERUN 
RETURN 

6 CALL ENDRW( NPLOT) 

END 

STOP 



114 



TITLE. THIRD ORDER POSITIONING SERVO 

* STEP RESPONSE 

* 

* MULTIPLE RUNS 

* 

* ANALOG AND DIGITAL SYSTEf-S 

* 

I NT EG RK3FX 
INTGER NPLOT, HOUR 

CONST NPL0T=4 , NCUR=1 , C=0 . 0 , D=1 . 0 , KV=6. 0 , RN=2. 0 , OP=l0 . 0 , Pl=0 . 2, . 

Z1=C.5,P2=12.5,Z2=1.25 

INITIAL 

KC=OP/ 6. 2831853 
TZ1=1.0/Z1 
TPl=1.0/Pl 
TZ2=1. 0/Z2 
TP2=1.0/P2 
DSRIVATIVS 

BIT=KC*C*D 

IG3IT=AINT(3IT) 

NT=RN-IG3IT-C*C 
VT=NT* ( D/KC+C ) 

VTN=VT*KV 

CD=LEDLAG ( 0 . 0 , TZ1 , TPl , VTN) 

COD=LEDLAG ( 0 . 0 , TZ2 , TF2 , CD) 

CDD=RS\LFL ( 0 . 0 , 0 . 2 , COD) 

CDOT=REALPL( 0.0, l.O.CDD) 

C=IN'TGRL ( 0 . 0 , C DOT ) 

CONTRL FINTIM— 20 . 0 , DELT=0 .0025, DSLS=0 .015 

PRINT 0.2,VT,CD0T,C,NT 

SAMPLE 

CALL DRV'G ( 1 , NCUR , TIME , C ) 

TERMINAL 

IF(NPLOT.EQ.l) GO TO 6 
NPLOT=NPLOT- 1 
NCUR=NCUR+1 

IF(KPLOT.E2.3) GO TO 2 
IF(NPLOT. EQ. 2) GO TO 3 
IF(NPLOT.B2.1) GO TO 4 

2 RN=9 . 0 
OP=45.0 
GO TO 5 

3 RN=36.0 
0P=180. 0 
GO TO 5 

4 RN=1. 25664 
Csl.o 
D=0.0 

5 CALL RERUN 
RETURN 

6 CALL ENDR/.'( NPLOT) 

END 
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